Quintom phase-space: beyond the exponential potential 
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We investigate the phase-space structure of the quintom dark energy paradigm in the framework 
of spatially flat and homogeneous universe. Considering arbitrary decoupled potentials, we find 
certain general conditions under which the the phantom dominated solution is late time attractor, 
generalizing previous results found for the case of exponential potential. Center Manifold Theory is 
employed to obtain sufficient conditions for the instability of de Sitter solution either with phantom 
or quintessence potential dominance. Another important issue is concerning the existence of a 
point corresponding to the standard quintessence dominated solution, which under certain condition 
on the potential, can mimic the dark matter behavior. This feature has important cosmological 
consequences to address de unified description of dark matter and dark energy in a single field. 
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I. INTRODUCTION 

Recent cosmological observations point to strong ev- 
idence for an spatially flat and accelerated expanding 
universe [H~0|- Despite the great agreement of observa- 
tions with the concordance model [4j] , it is a fact that 
quintom model whose Equation of State (EoS) can cross 
the cosmological constant barrier w = — 1 is mildly fa- 
vored 0, @| . A popular way to realize a viable quintom 
model and, at the same time, avoid the restrictions im- 
posed by the Non-Go Theorem 2 is the introduction of 
extra degrees of freedom. Following this recipe, the sim- 
ple quintom paradigm requires a canonical quintessence 
scalar field a and simultaneously a phantom scalar field 
<f> where the effective potential can be of arbitrary form, 
while the two components can be either coupled (Tl| or 
decoupled HO. 

The properties of the quintom models have been stud- 
ied from different points of view. Among them, the phase 
space studies, using the dynamical systems tools, are very 
useful in order to analyze the asymptotic behavior of the 
model. In quintom models this program have been car- 
ried out in [Iol - [l7j . In [lj[ the decoupled case between the 
canonical and phantom field with an exponential poten- 
tial is studied shown that the phantom-dominated scaling 
solution is the unique late-time attractor. In [ll[ the po- 
tential considers the interaction between the fields and 
shows that in the absence of interactions, the solution 
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1 The Cosmological Constant Model 

2 This theorem prohibits the crossing of the w = —1 boundary 
for models describe by perfect fluids or single scalar field with 
Lagrangian C= C{4>,dn4'd ,J '4 ) ) aI1 d coupled minimally to gravity 
in a FRW geometry [JujJ 



dominated by the phantom field should be the attractor 
of the system and the interaction does not affect its at- 
tractor behavior. This result is correct only in the case 
in which the existence of the phantom phase excludes 
the existence scaling attractors |l3[ . Some of this results 
was extended in [lj] for arbitrary potentials. In [ItJ the 
authors showed that all quintom models with nearly flat 
potentials converge to a single expression for EoS of dark 
energy, in addition, the necessary conditions for the de- 
termination of the direction of the w = — 1 crossing is 
found. 

The aim of this paper is to extend the study of Refs. 
[Tol [l2|-[l5j]-investigation of the dynamics of quintom 
cosmology- to include a wide variety of potential be- 
yond the exponential potential without interaction be- 
tween the fields, all of them can be constructed using the 
Bohm formalism [Tij|j2(jj of the quantum mechanics under 
the integral systems premise, which is known as quantum 
potential approach. This approach makes it possible to 
identify trajectories associated with the wave function of 
the universe [l8j when we choose the superpotential func- 
tion as the momenta associated to the coordinate field q. 
This investigation was undertaken within the framework 
of the minisuperspace approximation to quantum theory 
when we investigate the dynamics of only a finite number 
of models. Here we make use of the dynamical systems 
tools to obtain useful information about the asymptotic 
properties of the model. In order to be able to analyze 
self interactions potential beyond the exponential one we 
rely on the method introduced in Ref. [2l| in the context 
of quintessence models and generalized to several cos- 
mological contexts like: Randall-Sundrum II and DGP 
branes dHH, Scalar Field Dark Matter models 0, 
tachyon and phantom fields (26rl28| and loop quantum 
gravity H|. 

The plan of the paper is as follow: in section HT1 we in- 
troduce the quintom model for arbitrary potential and in 
section Mil we build the corresponding autonomous sys- 
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tern. The results of the study of the corresponding criti- 
cal points, their stability properties and physical discus- 
sion are shown in section Hv] Finally section |V] is devoted 
to conclusions. Finally, we include in the two appendices 
lAl and [B] the center manifold calculation of the solutions 
dominated by either the phantom or quintessence poten- 
tial (that is, the de Sitter solutions P5 or Pq). 



III. THE AUTONOMOUS SYSTEM 

In order to study the dynamical properties of the sys- 
tem ([2][5]) we introduce the following dimensionle ss p hase 
space variables to build an autonomous system [30, [HTj] : 



V6H' 



Xc/> = 



,y<7 = 



V3H 



(10) 



II. THE MODEL 



The starting action of our model, contai ning the canon- 
ical field a and the phantom field <j>, is [EL IllL H^: 



s = j d 4 xV=g Qi? - \gTd^Qu<T + V a (a)- 



+ -g^d^d^ + V^)), (1) 



A„ 



(11) 



Notice that the phase space variables \ a and A^ are 
sensitive of the kind of self interactions potential chosen 
for quintessence and phantom component, respectively 
and are introduced in order to be able to study arbitrary 
potential. Applying the above dimensionless variables 
to the system (J2KS]) we obtain the following autonomous 
system: 



where we used natural units (8irG =1) and V^ct) and 
V<p((()) are respectively the self interactions potential of 
the quintessence and phantom fields. From this action 
the Friedmann equations for a fiat geometry reads (Til 
El: 



2 J.2 



(3) 



and the evolution of the quintessence and phantom field 
are: 



dVJa) 

OCT 



(4) 



(5) 



where H = ^ is the Hubble parameter. Additionally me 
can introduce the total energy density and pressure as: 



PDE = Pa 



PDE = wpDE 



where 



w 



ct 2 -0 2 -K(ct)-T/ (0) 
ct 2 -ft + V a {a)+V^) 

and the dimensionless density parameter: 

r2/r , Qrh — 



PDE 



PDE 



and 



fl a + VLa, = 1 



(G) 



(7) 



(8) 



(9) 



^ = -3x a (1 + x% - xl) + f-ylK (12) 

dx<t, 
~dN 



-3x (l + x\ - x\) + 



+ -vl)K (is) 

~M = \ va i 6x2 ~ ^> x ° x ° ~ 6x i) ( i4 ) 

^ = -VEx a f(\ a ) (15) 
= _V6x 05 (A ) (16) 

where N = In a is the number of e-foldings and /(A CT ) = 
A 2 (r CT - 1) and c/(A ) = A 2 (r - 1) where: 



_ V a {a)dlV a {a) 



2 ' <t> 



(WW) 2 



(17) 



In order to get from the autonomous equation (I12j|16[) a 
closed system of ordinary differential equation we have 
assumed that the parameters T a and can be written 
as a function of the variables X a € K + and A^ € IR + 
respectively [2lj | . 

The phase space for the autonomous dynamical system 
driven by de evolutions of Eqs. (|12H16p can be defined as 
it follows: 



* = {{x a ,x^,y a ) :ya>0,xl-xl + yl< 1} 



x{(A ff ,A f )el 2 } (18) 



IV. CRITICAL POINTS AND STABILITY 

The critical points of the system (|12j|13|) are summa- 
rized in Table HI The eigenvalues of the corresponding 
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Jacobian matrices are show in Table |TTJ In both cases 
A* and arc the values which makes the functions 
f(K) = A 2 (r CT - 1) and g(X <t> ) = A 2 (T^ - 1) vanish re- 
spectively. In the same way we have chosen the notation: 



f'(K) 



dX„ 



dg(X cj> ) 



As we see from Table U the points P-f do not exist 
in the strict sense (x^ is purely imaginary at the fixed 
points). Point P5 is associated with a combination of 
a phantom potential whose first ^-derivative vanishes at 
some/several point/points, i.e., A^ = (this case include 
the exponential potential whose ^-derivative at any order 
vanish everywhere) and an arbitrary self interaction po- 
tential for the quintessence component (arbitrary value 
of A CT ). Point Pq is associated with a combination of a 
quintessence potential whose first er-derivative vanishes 
at some/several point/points, i.e., X a = (this case in- 
clude the exponential potential whose er-derivative at any 
order vanish everywhere) and an arbitrary self interaction 
potential for the phantom component (arbitrary value 
of A^). Point P7 is associated with a combination of a 
phantom potential whose first c6-derivative vanishes at 
some/several point/points, i.e., A^ = (this case include 
the exponential potential whose c6-derivative at any or- 
der vanish everywhere) and a self interaction potential 
for the quintessence component whose first c-derivative 
vanishes also at some/several point/points, i.e., A^ = 
. It is worth noticing that the existence of points P% 7 
P 3 ± , P4, P$ and P9 depends of the concrete form of the 
potential. From the table of the eigenvalues, notice, be- 
sides, that all the points belongs to nonhyperbolic sets of 
critical point with a least one null eigenvalue. 

Although all these critical points are shown in the Ta- 
bles U here we summarized their basic properties: 

• Pj 4 , P 3 ± , and P 2 ± correspond to a solution dom- 
inated by the kinetic energy of the scalar fields 
(stiff fluid solution: q = 2 and to = 1). The 
exact dynamical behavior differs for each points. 
P 1 corresponds to a phantom kinetic energy dom- 
inated (ft a- — and ftj, = 1). However, these 
points have a purely imaginary value of x^, thus, 
they do not exists in the strict sense. They 
have a three— dimensional center subspace and a 
two-dimensional unstable manifold (mi = 3 > 
0, ^.(1715) = 6 > 0). Thus they cannot be late-time 
attractors. P^ 1 represents an scaling regimen be- 
tween the kinetics energies of the quintessence and 
phantom fields (fi CT = l+arf and = —xV). These 
points depend of the form of the potentials and un- 
der certain conditions they have a four dimensional 
unstable subspace which could correspond to the 



past attractor. However, this point is unphysical 
since Q$ < . P^ is dominated by the quintessence 
kinetic term (f2 CT = 1 and f2^ = 0). Since they 
are non-hyperbolic due to the existence of two null 
eigenvalues, we are not able to extract information 
about their stability by using the standard tools 
of the linear dynamical analysis. However, since 
these points seems to be particular cases of P 3 ± , 
they should share the same dynamical behavior. 
Because all of these points are nonhyperbolic, as we 
notice before, we cannot rely on the standard lin- 
ear dynamical systems analysis for deducing their 
stability. Thus, we need to rely our analysis on nu- 
merical inspection of the phase portrait for specific 
potentials or use more sophisticated techniques like 
Center Manifold theory. 

• P4 is an scaling solution between the kinetic and 
the potential energy of the quintessence component 
of dark energy. This solution in sensitive to the 
explicit form of the potential. This is always a 
saddle equilibrium point in the phase space since 
m 2 — (A*) 2 and m 4 = i((A*) 2 — 6) are of oppo- 
site sign in the existence region of this point. It 
represents an accelerated solution for a potential 
V a (a) whose function f(\ a ) vanish for X a = A* 
in the interval — \[2 < A* < a/2, leading to a 
— 1 < w e // < —1/3. When A* = the critical point 
P4 becomes in P 6 . In the regions —a/6 < A* < — v2 
or \/2 < A* < a/6, the critical point P4 represents 
a non-accelerated phase. A very interesting issue 
of this critical point appears when, for an specific 
form of the quintessence potential, A* = ±vo, driv- 
ing to w e ff = 0. This means that the quintessence 
field is able to mimic the dark matter behavior and 
this could be a nice scenario to address the unified 
description of the dark matter and of the dark en- 
ergy in a single field. This result is very important 
in order to describe the dark matter evolution by 
mean of an specific quintessence scalar field that 
fulfill the above condition due to the matter (dark) 
domination era is a transient stage of the cosmic 
evolution. 

• P5, Pe and P7 represents solutions dominated by 
the potential energies of the potentials (all of them, 
but P7, represent de Sitter solutions: q — —1 and 
w e ff = — 1). Once again the exact dynamical na- 
ture differs from one point to the other: P5 is domi- 
nated by the potential energy of the phantom com- 
ponent (0 CT = and fl^ — 1). Because of the exis- 
tence of two null eigenvalues is not possible to con- 
clude about its dynamics. However it has a three- 
dimensional stable manifold for g(0) < (in the 
interval g(0) < — | it has to complex conjugated 
eigenvalues with negative real parts). In this cases 
it is worthy to analyze its stability using the center 
manifold theory. P§ is a critical point dominated by 
the quintessence potential energy term (fi CT = 1 and 
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TABLE I: Properties of the critical points for the autonomous system (|12II16I 



Label 




y<r 




Act 


X<f, 


Existence 


Oct 






q 


Weff 


Pt 








±i 


Act 




Non real 





1 




2 


1 


Pt 


±i 








Act 


X<l> 


Always 


1 







2 


1 


Pt ±, 









A* 


Kj, 


11 


1 + xl 


-4 




2 


1 


F 4 




V 6 





A* 




-V6 < A* < V& 


1 







+ (A i )2 

T 2 


i+ (A i )2 

' 3 


P 5 











Act 





Always 





1 




-1 


-1 


Pe 





1 








A^ 


?j 


1 







-1 


-1 


Pt 





2/o- 











< y a < 1 


vl 




2 


-3y CT 


1-22/ct 


Ps 








a ; 

V6 


Act 


A^ 


Act G K 





i 


-1 


(AJ) 2 
2 


1 3 



TABLE II: Eigenvalues of the critical points 



Label mi m.2 tn3 



pf 3 o o t»\/6<?'(a;) 6t»\/6a; 

6 tV6/'(Act) 3+^/1^ 



Pt -V6c/(A;)^ =F>/6/'(A;)^ + 1 3 T ^/|^ + 1Act &-^x^X% 

P 4 —f'{\%)\a* (\a*f ±((\a*f-6) § ((Aa*) 2 - 6) 



P 5 -3 -| (^l + ^l + | ff (0)j -| (^1-^1 +§,,(0) 

P 6 -3 ^1 - |/(0)) _|(l_^l_|/(0)) 

Pt \ (-^9 - 12/(0)2/1 - 3) 1 (^9 - 12/(0)y 2 - 3) § (-^9 - 12 S (0) (y 2 - 1) - 3) | (^9 - 12 fl (0) (^ - 1) - 3) 

A <?'(a;)a; -H>4) 2 -|(( A ;) 2 + e) -I ((a;) 2 + e) 



Clf = 0), despite its nonhyperbolicity, it has three- 
dimensional stable manifold for /(0) > (in the 
case /(0) > I it has to complex conjugated eigen- 
values with negative real parts), thus, it is wor- 
thy to analyze its stability using the center mani- 
fold theory. P7 denotes a segment (curve) of non- 
isolated fixed points, representing a scaling regimen 
between the quintessence and phantom potential 
(Q a = y 2 a and = 1 — y^). The existence of one 
non-zero eigenvalue is due to the fact that it is a 
curve of fixed points. As an invariant set of non- 
isolated singular points it is normally-hyperbolic, 
since the eigenvector associated to the zero eigen- 
value, (0, 0, 1, 0, 0) T , is tangent to the curve. Thus 
its stability is determined by the sign of the re- 
maining non-null eigenvalues. Hence, it is stable 
for < y a < l & /(0) > & g(0) < or a saddle 
otherwise. 



• Ps is a line of fixed points parameterized by X a € M. 
The existence of one non-zero eigenvalue is due to 
the fact that it is a curve of fixed points. As an 
invariant set of non-isolated singular points it is 
normally-hyperbolic, since the eigenvector associ- 
ated to the zero eigenvalue, (0, 0, 0, 1, 0) T , is tan- 
gent to the curve. Thus its stability is determined 
by the sign of the remaining non-null eigenvalues. 
From table Hi] follows that Ps admits a four dimen- 
sional stable subspace provided g'(ApA^ < 0, thus, 
the invariant curve is stable. It represents acceler- 
ated solutions dominated by the phantom potential 
providing a crossing through the phantom divide 
= and = 1). For every value of A J this 
point provide the typical superaccelerated expan- 

(A* ) 2 

sion of quintom paradigm (w = —1 — — | — ) the 
only exception occurs when AJ = recovering the 
behavior of the de Sitter solution P 5 (ui = — 1). 



5 



This line of critical point corresponds to the stable 
point P in [l2| and B in (To| (phantom dominated 
solution). In our case the line P% is the late time 
stable attractor provided g'(XV)X% < 0, otherwise, 
it is a saddle point. 

Finally, in order to examine the stability of the nonhyper- 
bolic points that cannot consistently be studied via the 
present linear analysis, we present a concrete example. 
We provide a numerical elaboration of the phase space 
orbits of the corresponding quintom model. 



A. V(a, </>) = Vo sinh 2 (aa) + Vi cosh 2 (^) 

This potential is derived, in a Friedmann-Robcrtson- 
Walkcr cosmological model, from canonical quantum cos- 
mology under determined conditions in the evolution of 
our universe 3 , using the bohmian formalism fl8| . For this 
potential: 

/(A CT ) = -^+2a 2 , \l = ±2a, /'(A*) = —A* (19) 



and 



\ 2 



9(*4,) = —* + W, a; = ±2/3, g'(\;) = -x; (20) 

from Table |TT] and equation (I2TJ1) we see that condition to 
ensure that Point Ps has a four dimensional stable sub- 
space is always satisfied due to the opposite sign between 
AJ and g'(X%). In order to having achieved success scalar 
field dark matter domination era we need that A* = ±^/3 
so the potential parameter a can be adjusted using (TK 
We also impose one of the following conditions: 




FIG. 1: Trajectories in phase space (x a , y&, y^) for dif- 
ferent sets of initial conditions for potential V(a, tf>) = 
Vo sinh 2 (out) + Vi cosh 2 (/30). The free parameter have been 
chosen to be (a, 0): (— 1/3/2, 0.35). This parameter selec- 
tion guarantee that point P4, black point in this graphic, rep- 
resents an scalar field matter dominated era with a typical 
saddle dynamics. Late time attractor is the phantom field 
dominates solution (Ps in Table [TTJ) . 




or 



a* = \/3 & Ai < -Ve & i< x a < 



\ 



G 



FIG. 2: Trajectories in phase space (x a , y a with the same 
parameter selection of Fig. \T\ Critical point P4 is represented 
by black point in this graphic, corresponding to an scalar field 
matter dominated solutions 



V3. 



\/6<A:<0&l<x CT <%/2 



CONCLUSIONS 



or 



A! 



-V3k A; < -V6k 



\ 



< x a < -1 



-V3 & - VE < XI < & - V2 < x„ < -1. 



to guarantee that the stiff matter type solution of the 
quintom cosmology be the past-attractor. 



This is part of a forthcoming paper. 



In the present paper a thorough study of the phase 
space of quintom model has been undertaken. The results 
are valid for those potential, without interaction between 

the fields, for which the quantities r CT 



VA<T)diVa(<y) 



and 



I = (q^v ((fi))' 1 can written as a function of the 
variables A CT and A^. 

It has been found that for g'(X^)X*^ < 0, the late 
time attractor are always the phantom dominated so- 
lution (P s ) generalizing the result shown in [l(| [3 for 
exponential potentials (but more generally, for potential 
satisfying X a w const and A^ ~ const) . Otherwise, it is a 
saddle point. The Universe evolves from a quintessence 
dominated phase to a phantom dominated phase crossing 
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FIG. 3: Trajectories in phase space (y„, y$) with the same 
parameter selection of Fig. [T] The accelerated de Sitter solu- 
tion Pt, dashed line, is a transient era in the evolution of the 
Universe being the late time attractor the phantom field dom- 
inated solution, black point in this figure, allowing a crossing 
through the phantom divide. 



the w e ff = — 1 divide line as a transient stage 32 1. 

Center Manifold Theory have been employed to ana- 
lyze the stability of dc Sitter solution either with phan- 
tom (P5) or quintessence potential dominance (Pe)- Af- 
ter deriving the evolution equation on the center mani- 
folds and making several numerical integrations we have 
concluded that in both cases the corresponding de Sitter 
solution is unstable (saddle-like). For P 5 we have used an 
analytical argument, whereas for Pq our conclusion was 
supported on numerical experimentation. 

Another important issue is concerning the existence of 
a point, P4, corresponding to the standard quintessence 
dominated solution, which under certain condition on the 
potential, can mimic the dark matter behavior. This fea- 
ture has important cosmological consequences to address 
de unified description of dark matter and dark energy in a 
single field. The stable character of P4 have been clearly 
illustrated by resorting to phase plane diagrams for the 
potential obtained from a canonical quantum cosmology. 



Appendix A: Center manifold dynamics for the 
solution dominated by the potential energy of the 
phantom component Pg 

In this section we will shown how we can apply the 
center manifold theorem to study the stability of non- 
hyperbolic point P5 corresponding to the solution dom- 
inated by the potential energy of the phantom compo- 
nent [33l ] . First, we restrict our attention to the domain 
~| < 3(0) < to dealing with real eigenvalues. The 
first step is to translate the point P5 (x a = 0, x^ = 
ijcr = 0, A CT = a, = 0) to the origin, where a denotes 
an arbitrary value for \ a . The next step is to transform 



the system to its real Jordan form: 

u = Zu + P(u,v) 
v = Pv + G(u,v) 



(Al) 
(A2) 



where the square matrices Z, P have 2 zero eigenvalues 
and 3 eigenvalues with negative real part, respectively. 
In order to do that we introduce the new variables: 



u 2 = -\j-x a f{a) - a + A 



Vl 



V 2 = 



3^/(0), 

2V6g(0)x^ + (Vl2 5 (0) + 9 - 3) A 



2^12*7(0) + 9 



v/12 ff (0) + 9 + 3 A - 2\/6 5 (0)^ 

v 3 = > (A3) 

2Vl2<?(0) + 9 



Using the above transformation, the system (| AIM A2|) is 
given explicitly by: 



Ui = F 1 (ux,U 2 ,V 1 ,V2,V 3 ) 

3«i(/(a + u 2 + ui) - /(a)) 



U 2 = 



/(a) 

EE F 2 {ui,U 2 ,Vl,V 2 ,V 3 ) 

vx = -3ui + G 1 (u 1 ,u 2 ,v 1 ,v 2 ,v 3 ) 
v 2 = - ^-^12^(0) + 9 - 3) v 2 + G 2 (u 1 ,u 2 ,v 1 ,v 2 ,v 3 ) 



(A4) 

+ H(ui,u 2 ,v 1 ,v 2 ,v 3 ) 

(A5) 
(A6) 



(A7) 

v 3 = - (V 12 3(°) + 9 - 3) v 3 + G 3 (u 1 ,u 2 ,v 1 ,v 2 ,v 3 ) 

(A8) 

near the non- hyperbolic fixed point P5 where Pi , H, G± 
... G 3 are homogeneous polynomials in the coordinates 
(ui,u 2 ,vi,v 2 ,v 3 ) of degree greater than 2. Following the 
standard formalism of the center manifold theory, the 
coordinates which correspond to the non-zero eigenvalues 
(«l, i>2, ) can be approximated by the functions: 

ki(ui,u 2 ) — a\u\ + a 2 u\ + a 3 uiu 2 + a 4 ulu 2 + a 5 ul+ 
+ a e u x ul + a 7 ul + ... + 0{u1,u 2 ) (A9) 

k 2 (ui,u 2 ) = biuj + 6 2 m'i + b 3 uiu 2 + b^u\u 2 + b 5 ul + 

+ b 6 u 1 ul + b 7 ul + ... + 0(u^,u^) (A10) 

k 3 (ui,u 2 ) = ciu\ + c 2 u\ + c 3 uiu 2 + c A u\u 2 + c^u\+ 

+ C 6 Ui«2 + C 7 U% + ... + (9(u™,U2) (All) 

with this set of functions we can solve, to any n desired 
degree of accurancy, the quasilinear partial differential 
equation for the center manifold: 

Dk(u) [Zu + P(u, fc(u))] - Pfc(u) - G(u, fc(u)) = 

(A12) 
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In our case: Z = 



P 








and 



/-3 

i(-yi2^0j" 

y 



9-3 



t — > +oo, unless a 
unstable. 



0. Thus, for -f < 0(0) < 0, P 5 is 



- there are two complex eigenvalues. In 



For 5 (0) < 4 

\ this case, in order to obtain the real Jordan Form, we 



(v/12.g(0) + 9-3) j 



introduce the new variables 

V 2 + V 3 



, v 3 



V2 ~ V3 

2i ' 



fc(u) = 



' ki(ui,u 2 ) 
k 2 (ui,u 2 ) 
\ k 3 (ui,u 2 ) 



Using the above transformation, the system (| AIM A2|) is 
given explicitly by: 



Dk(u) 



I dki dk± 

du\ du2 
dki 
du\ du2 
dks dk 3 
du\ du2 



G(u,fc(u)) 



1 G\ (Ui,U 2 ,ki(Ui,U2),k2(u-i,U2),k 3 (ui,U2)) 

G 2 (ui,u 2 ,k 1 (ui,u 2 ),k 2 (u 1 ,u 2 ) 1 k 3 (u 1 ,u 2 )) 
\ G 3 (ui,U2,k 1 (ui,u 2 ),k 2 {u 1 ,u 2 ) 1 k 3 (u 1 ,u 2 )) 



ui 
u 2 ' 

V 2 ' 



Fl(ui,U2,Vl,V2,V 3 ) 

F 2 (ui,u 2 ,vi,V 2 ,V 3 ) 

-3vi +Gi(u 1 ,u 2 ,vi, V 2 ,V 3 ) 
1 



(A19) 
(A20) 
(A21) 



--Vt - - V-12<?(0) - 9V 3 + G 2 {ux, m, «i, V 2 , V 3 ) 

(A22) 

-1^-123(0) - W 2 - ~V 3 + G 3 {ux,U2,v 1 ,Vi,V 3 ) 

(A23) 



In order to solve equation (|A12[) we put together Z, P, 
k(ui,u 2 ), Dk(u\,u 2 ) and G(ui,u 2 ,k(ui,u 2 )), then we 
equate equal powers of u\ and u 2 , and in that way we 
compute k{u\,u 2 ). Finally we obtain, in the neighbor- 
hood of P5 , the reduced system: 



u = Zu + F(u,k(u)). 



(A13) 



The above procedure applied to equations (|A4HA8[) leads 
to: 

vi = \af{a)u\ + i/(o)«?«a + 0(4); v 2 =v 3 = 0(4). 

(A14) 

Neglecting the fourth order terms, the evolution equa- 
tions on the center manifold are 



1 



2„,3 



-a u 



-u\{a + u 2 )f(a) 



(A15) 
(A16) 



The orbit of (|A15|) - (|A16[) passing through (uio,u 20 ) is 
given by 



"2 



,T>0, 



1 



a if f(a) ^ ^ ^ Al7 ^ 
if/(o) = ' 

(A18) 



For f(a) > 0, the orbits approach the point with co- 
ordinates (ui =0,u 2 = —a) when r — » +00. If /(a) < 0, 
then, as r — > +00, u\ tends to zero and u 2 becomes un- 
bounded. Generically, the origin is not approached as 
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where F\,G\ ... G 3 are homogeneous real polynomials in 
the coordinates (ui,u 2 , v\, V2, V 3 ) of degree greater than 
2. Usig the same procedure as before we obtain that the 
center manifold is given locally by the graph 

vi = \af{a)u\ + ^f{a)u\u 2 + 0(4); V 2 = V 3 = 0(4). 

(A24) 

Thus the dynamics on the center manifold is given by 
the system (|A15l) - (|A16p analyzed before. Sumarizing, 
for 3(0) < 0, P5 is unstable. 



Appendix B: Center manifold dynamics for the 
solution dominated by the potential energy of the 
quintessence component Pe 

In this section we will shown how we can apply the 
center manifold theorem to study the stability of non- 
hyperbolic point Pq corresponding to the solution dom- 
inated by the potential energy of the quintessence com- 
ponent [33j. 

The first step is to translate the point P 6 (x a = 0, 
%cf> = Do = 1, A CT = 0, = a) to the origin, where a 
denotes an arbitrary value for A^. The next step is to 
transform the system to its real Jordan form: 



u 

V 



Zu + P(u,v) 
Pv + G(u,v) 



(Bl) 
(B2) 



where the square matrices Z, P have 2 zero eigenvalues 
and 3 eigenvalues with negative real part, respectively. 
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In order to do that we introduce the new variables: 
ui = -2a(y a - l)g(a), 

u 2 = -g5(o) (y6x<p — 2a{y a — 1)J - a + \<p, 
vi = -g(a) (ySxcf, - 2a(y a - 1)) , 

2V6.f(0)x a + U$ - 12/(0) - 3) A ff 

v 2 = — , 

2^/9-12/(0) 
-2^/(0)^ + (v/9- 12/(0) + 3) A fl 



v 3 



2^/9-12/(0) 



(B3) 



Using the above transformation, the system (|BlHB2j) is 
given explicitly by: 

Ux = Fi(ui,U2,V 1 ,V2,V 3 ) (B4) 

, 3ui(5'(a + 'U2 +v±) — g(a)) uig(a + u 2 + vi) 

u 2 = + 

g{a) g{a) 

+ H(ui,u 2 ,Vi,V2,v 3 ) ee F 2 {ui,u 2 ,vi,v 2 ,v 3 ) (B5) 
vx = —3vt + Gi(ut,u 2 ,vi,V2,Va) (B6) 

v 2 '= i (-v/9 -12/(0) -3) wa + GaCux.ua.wi.wa.ws) 

(B7) 

«a' = i (v/9 -12/(0) -3) «3 + G 3 («i,tia,t;i,«2,«3) 

(B8) 

near the no n- hyperbolic hxed point Pg where Fi,H,Gi 
... G 3 are homogeneous polynomials in the coordinates 
(ui,u 2 ,vi,v 2 ,v 3 ) of degree greater than 2. Following the 
standard formalism of the center manifold theory, we ob- 
tain that the center manifold of P§ is given by the graph 



Vl - TE^oj + 12a S (a) +9^ ~ + C7 ^- ) ' 

«2=0(3),«3 = 0(3). 



(B9) 



Neglecting the third order terms, the evolution equations 
on the center manifold are 



(B10) 
(Bll) 



.9(a) 

Ul u 2 (2a 2 + 3) u\ 
a Y2ag{a) 



Let us assume g(a) / 0,a / 0. Hence, the orbit of 
(IB10|) - (|B11|) passing through (wio,U2o) is given by 



1*2 = 



((2a 2 + 3) u w + 12a 2 u 2 o - ^2u 2 og(a)) ( u\ 



12 (a 2 - 5 (a)) 



"10 



(2a 2 + 3) mi 
+ 12 (.g(a)- a 2 )' 
ui g(a) 



u w g{a) - au W T 



T > 0, 



(B12) 
(B13) 



In order to investigate the stability of the center manifold 
of Pq we have resorted to several numerical integrations of 
the system (|B10j) - (|Bllj) . We find four typical situations 
that suggest that Pq is unstable (saddle type). 




Vector field in the plane 

2/ \ , T/ i2| 



U2) for the potential 



FIG. 

V{<7, (j)) — Vo s'mh' 1 (aa) + Vi cosh^(/30). The free parameter 
have been chosen to be (a, /?, a): (— \/3/2, 0.35, 0.30). In 
this case g(a) = 0.20 > 0. The sign of u\ is invariant. For 
u\ < the origin is approached as the time goes forward 
whereas for u\ > the orbits departs from the origin. Thus 
The accelerated de Sitter solution Pa is a transient era in the 
evolution of the Universe. 

In the figure 0] is displayed the vector field in the 
plane (u\, u 2 ) for the potential V(cr, <f>) = Vq sinh 2 (a<r) + 
Vi cosh 2 (f34>) ■ The free parameter have been chosen 
to be (a, ft, a): (— s/3/2, 0.35, 0.30). In this case 
g(a) = 0.20 > 0. The sign of u% is invariant. For 
Ui < the origin is approached as the time goes for- 
ward whereas for u\ > the orbits departs form the 
origin. For the choice (a, (3, a): (—-s/3/2, 0.35, -0.30). 
we have g(a) = 0.20 > 0. The figure is similar to |4] with 
the arrows in reverse orientation. Thus, this numerical 
elaboration suggest that the accelerated de Sitter solu- 
tion Pq, (the origin of coordinates), is a transient era in 
the evolution of the Universe for g(a) > irrespectively 
the sign of a. 

In the figure [5] is represented the vector field in the 
plane (ui, u 2 ) for the potential V(a, tf>) = Vo sinh 2 (ever) + 
Vi cosh 2 (/3(/>). The free parameter have been chosen 
to be (a, P, a): (— v3/2, 0.35, 1.30). In this case 
g(a) = —0.60 < 0. The sign of u\ is invariant. All 
the orbits departs from the origin. Thus the accelerated 
de Sitter solution Pq is a transient era in the evolution 
of the Universe. For the choice (a, (3, a): (—■ y/3/2, 0.35, 
—0.30). we have g(a) = 0.20 > 0. The figure is sim- 
ilar to H] with the arrows in reverse orientation. Thus, 
this numerical elaboration suggest that the accelerated 
de Sitter solution Pq, is a transient era in the evolution 
of the Universe for g(a) < irrespectively the sign of a. 

As in the appendix^ for analyzing the case of complex 
eigenvalues, we can introduce the new variables 



V 2 



v 2 + v 3 



, v 3 



v 2 - v 3 
2i 
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FIG. 5: Vector field in the plane (ui, U2) for the potential 
V(a,(f>) = Vo sinh 2 (aa) + Vi cosh 2 (/30). The free parameter 
have been chosen to be (a, /?, a): (— V3/2, 0.35, 1.30). In this 
case g(a) — —0.60 < 0. The origin is of saddle type. Thus 
the accelerated de Sitter solution Pe is a transient era in the 
evolution of the Universe. 



for deriving the real Jordan form of the Jacobian. The 
procedure is straightforward, so we won't enter into the 
details here. 
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